Wavefront solution in extended quantum circuits with charge discreteness 
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A wavefront solution for quantum (capacitively coupled) transmission lines with charge discrete- 
ness (PRB 64, 235309 (2001)) is proposed for the first time. The nonlinearity of the system becomes 
deeply related to charge discreteness. The wavefront velocity is found to depend on a step disconti- 
nuity on the (pseudo) flux variable, /, displaying allowed and forbidden regions (gaps), as a function 
of /. A preliminary study of the stability of the solutions is presented. The dual transmission line 
hamiltonian is proposed and finally, we find a connection with the (quantum) Toda lattice. 
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I. INTRODUCTION 

Today, the broad field of nanostructure is at the heart 
of many technological deviceaiSiSiii^. Naturally, at this 
scale and for low temperature, quantum mechanics plays 
a fundamental role ^. Recently, much effort has been 
devoted to study nanostructures, using as model that of 
quantum circuits with charge discretenessSiStSiifliiLi^. 
Work has been published on subjects such as persistent 
currents coupled quantum circuita^iS, the electronic 
resonatorifl.^ quantum point contactsii, and others-. In 
this article, we are interested in spatially extended quan- 
tum circuits (transmission lines) with charge discrete- 
ness. This is an interesting specific theoretical subject of 
broad potential applications, since nano-devices could be 
put together forming chains, and then be viewed as elec- 
tric transmission lines. For instance, electric transport 
properties in DNA have been measured recentljiiSii^ii^; 
some degree of disagreement related to conducting prop- 
erties exists, nevertheless, it is clear that the DNA 
molecule could be viewed and modeled as a quantum 
transmission line. Moreover, molecular electronic cir- 
cuits'^ are actively studied theoretically and experimen- 
tally. In such systems, chains of individual molecules 
form a line of circuits. Therefore, quantum circuits is 
a broad field involving future applications from the per- 
spective of nano-devices, electric transmission in macro- 
molecules, left-handed materials, Toda lattice, and oth- 
ers. 

In this work we will consider a wavefront solution for 
an extended quantum circuit (transmission line). For 
the first time, a sequence of bands and gaps is found 
and characterized for this specific system with charge dis- 
creteness, giving hopes for the description of more com- 
plex extended systems (dual transmission lines, complex 
basis, etc.). 

In section^]we will introduce the Hamiltonian for cou- 
pled circuits, and the equations of motion for the spa- 
tially continuous systems. In Sec. IIIII the wavefront 



solution is obtained and the band-gap structure is char- 
acterized, in Sec. llVl the stability question is considered. 
In Sec. the dual line hamiltonian is introduced. In 
Sec. I VII the connection with Toda lattice is presented. 
Finally, we state our conclusions. 



II. QUANTUM TRANSMISSION LINES WITH 
DISCRETE CHARGE 



It is known that for a chain of quantum capacitively- 
coupled quantum circuits with charge discreteness ((je), 
the Hamiltonian may be written as^: 



H = 



1 — — OC 



2r 



(1) 



where the index I describes the cell (circuit) at position 
I, containing an inductance L and capacitance C. The 
conjugate operators, charge Q and pseudoflux 0, satisfy 
the usual commutation rule [Qi, (jji'] = ihSi^i' . A spatially 
extended solution of Eq. will be called a cirquiton like 
solution, corresponding to the quantization of the classi- 
cal electric transmission line with discrete charge (i.e. el- 
ementary charge Qe). Note that in the formal limit 
the above Hamiltonian gives the well-known dynamics re- 
lated to the one-band quantum transmission line, similar 
to the phonon case. The system described by Eq. |^ is 
very cumbersome since the equations of motion for the 
operators are highly nonlinear due to charge discreteness. 
However, this system is invariant under the transforma- 
tion Qi — > (Qi + a), that is, the total pseudoflux operator 
^ — commutes with the Hamiltonian; in turn, the 

use of this symmetry helps us in simplifying the study of 
this system. 

To handle the above Hamiltonian, we will assume a 
continuous approximation (infrared limit); that is to say, 
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we shall use (Q/ — Qi-i) ~ dQ/dx^-. In this way, it is 
possible to re-write the sum of Eq. as the integral 

oc oc 

^ « J dx, where a; is a dimensionless variable used 

/— — oo — oo 

to denote the position in the chain (i.e. the equivalent of 
I), since we have not introduced so far the cell size of the 
circuit. In this approximation, the Hamiltonian becomes 



of solutions: a projection operator (/) — P, satisfying 
= P, and (p = a satisfying = 1. For simplicity, 
we shall consider only the first case (a projector). Con- 
sider the pseudoflux operator only in one LC-cell of the 
chain, and its spectral decomposition in the Schrodinger 
picture (pceii = J >< 4'\ d(f>- Now, pick-up only one 
term from there (call it 00, say) and consider now the 
well-defined operator 



H= Hdx = 



2h 



dQ 



dx, 



(2) 



where Ti represents the Hamiltonian density operator 
for the fields 0(a;) and Q{x) and where [Q{x),(t>{x')] = 
i^Sx,x'- From the above Hamiltonian we find the equa- 
tions of motion (Heisenberg equations): 



d - 

dt 



- — sin( — ^ 
Lqe h 



(3) 
(4) 



III. WAVEFRONT SOLUTIONS 



As stated previously, we are interested in wavefront- 
like solutions of the system of Eas. H3l4|l . We proceed 
in the standard way^, by assuming travelling wave so- 
lutions for our operators, i.e., we define a new variable 
z = X — vt, and assume 



Q{x,t) = Q{z), 
(t>{x,t) = 4){z), 



(5) 
(6) 



where v stands for the unknown propagation velocity 
with dimensions of inverse time, since we did not intro- 
duce a grid length parameter in our continuum approxi- 
mation. Therefore, from the Heisenberg equations of mo- 
tion 1)314(1 . the wavefront equations in the new variable z 
become 



dz 
dz 



— sin — ( 

Lqe h 



(7) 
(8) 



From the above pair of equations we obtain a closed 
equation for the pseudoflux operator resulting in the 
"eigenvalue" problem: 

^ • '^^{z)=v'^{z), (9) 



LCqe h 



where the integration constant has been chosen as zero 
by simplicity (however, see eg. ITU) . Equation ^ corre- 
sponds to an eigenvalue problem for the non-linear super- 
operator C{(j)) = sin(0). There is at least two kinds 



= fPo where Pq = |0o >< 



(10) 



where / is an arbitrary pseudoflux parameter, and re- 
place into Eq.( O. At this point, a note concerning 
the validity for the commutation rules may be in order: 
the projector operator Pq defines a one dimensional sub- 
space in which subspace the commutation rule between 
(j) = f Po and the projected charge operator PqQPq, is 
verified. 

Now, since Pq is a projector, then the equation for the 
pseudoflux / becomes related to the velocity by 



1 sin{qef/h) 
LC (qef/h) ■ 



(11) 



Since both signs (±/) are possible, then we can construct 
the wavefront solution (step 2/) of the equations of mo- 
tion (EEj): 



4>sol{z) 



-fPo, z<0 



(12) 



corresponding to a solution with zero total flux (Sec. QJ. 
Concerning the matching condition at z = 0, the solution 
()12|l is in complete agreement with the matching implying 
exphcitly Eas. lf7IS|l . 

The condition > on the wave-front velocity gives 
the band-gap conditions on the system. In fact, from 
(|ll|l the restriction 



9e f 

h J 



> 0, 



(13) 



means that there exists a sequence of bands and gaps. 
The figure Q shows a plot of the wavefront velocity 111|) 
for different values of the pseudoflux /, in which it is seen 
an alternating sequence of bands and gaps, corresponding 
to propagating and forbidden modes. Note that the main 
allowed band {—wh/qe < / < Trh/qe) is twice as wide as 
the other allowed bands. 

Recall that the integration constant was set equal to 
zero to obtain @, now consider the nonzero case, then 
the wavefront velocity becomes formally 



LCqe 



— sm — G 
n 



(14) 



where C is an integration constant. 
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FIG. 1: Plot of the velocity of the wavefront, as a function 
of the flux parameter /. As specified by Eg.f IllH there is a 
structure of bands and gaps. The main allowed velocity band 
is twice as wide as the other allowed bands. This structure 
is a direct consequence of charge discreteness, an effect that 
dissapears in the limit —> 0. 



IV. STABILITY 



Since = fPo, and as said Pq is a projector, then from 
the above equation (and (3-4)) we have the linear evolu- 
tion equation for the perturbation e: 



92 ^ 92 



£ e 



(17) 



We shall consider two cases: 



(a) e = e(x,t)Po- Here [e, Pq] = 
linear wave equation LC-^- 



the perturbations arc unstable when (cos ^ 



and then we have the 
-e — {cos -^e and 

' < 0. 



(b) £ = e{x,t){I 
becomes LC 
stable. 



Pq). In this case the wave equation 
— ■§^s, and the perturbation is 



9*2 = 



We present some preliminary results concerning the 
stability of the solution UlUI) ; a more detailed treatment 
is currently under consideration. Consider now the per- 
turbative solutions 



(f) = (l>sol + £ 

Q = Qsol + rji 



(15) 



where the operators e and 77 are the perturbation and 
then, assumed with small eigenvalues. Moreover, con- 
sider the well known perturbative expansion (— ) 



Re 



,i<t> 



e e 



(16) 



V. DUAL TRANSMISSION LINE 

It is well known that the direct classical L-C transmis- 
sion line (related to equation ^ has associated a dual 
transmission line: in the direct line the interaction be- 
tween cells is through capacitances, while in the dual it 
is through inductances. The dual line is closely related to 
the so-called left-handed materials, and then its quantiza- 
tion is actually important. Moreover, the role of charge 
discreteness must be also considered. This two step pro- 
cess (quantization and charge discreetness) could be real- 
ized in analogy with the direct line (section P) but in this 
case long range interactions between cells appear in the 
Hamiltonian. In fact, the expression for the Hamiltonian 
is (!«): 



H = 



2nLql 



■E 



^ik{l-n) 

dk sin — ( 

1 — cos k 2n 



>i sm 



2h' 



2C 



(18) 



The equation of motion the charge and pseudoflux may 
be obtained after some algebra. We shall not continue the 
description of the dual line, but we shall only mention 
here that the expressions are quite involved and, so far, 
no explicit solution have been found. 



VI. TODA LATTICE AND CIRQUITONS 

The cirquiton Hamiltonian |^ gives the equations of 
motion for charge and pseudo-flux in the cell I : 



d 1 

'I = — { Qi+i + Qi-i - "^Qi 



dt 



c 



(19) 



4 



dt 



Qi 



— sm— 0i. 



(20) 



Equations H3I4() in section ^ are the spatially continuous 
version of the above pair of equations. Then, for the 
pseudo-flux operator we have the closed non linear equa- 
tion of motion: 



LCt7ttt0; = — sm— +sm— - 2sm— 0z 
\ h n n 



(21) 

which, in the formal limit Qe —>■ gives the usual 
one-band system with frequency spectrum uj{k) = 
sin(fc/2)|, as expected. Consider now the equation: 



/LC 



+ e- 



2e- 



(22) 



then any hermitian solution of this last equation is also 
solution of the cirquiton equation (|21|l . To show this, it 
suffices to conjugate Eq. (I22|) and subtracting. Notice 
that any hermitian solution of Eq. (|22|) is also a solution 
of l|21() . but the inverse is not necessarily true. 

On the other hand, the broad field related to the 
Toda lattice is studied in depth in a variety of branches 
in physicsii including non-linear physics, statistical me- 
chanics, classical electric circuits, etc. The classical evo- 
lution equation for the Toda lattice has the generic form: 



2e 



-B4>i 



) (23) 



where M,A and B are real constant. So, we have the 
important result: the formal replacement qe — > iqe trans- 
forms 122|l into the (quantum ) Toda lattice equation. 
This makes a direct connection between Toda lattice and 
cirquiton theory. 



VII. 



FINAL REMARKS 



For the quantum electric transmission line with charge 
discreteness described by the Hamiltonian Eq.(^, and 
equations of motion (|3I4|I . a one parameter (/) wave- 
front solution was found (Eqs. I11I12|I . The condition 
Ea. Hll|l on the velocity generates a band-gap structure 
dependent on the pseudoflux parameter / (see figure 
namely, there exist regions (values of /) for which a soli- 
tary wavefront propagates with constant speed. The ex- 
istence of the band-gap structure described is the main 
result of this work. 
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